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The spherical symmetry Black holes are considered in expanding background. The singularity 
line and the marginally trapped tube surface behavior are discussed. In particular, we address the 
conditions whether dynamical horizon forms for these cosmological black holes. We also discuss 
about the cosmological constant effect on these black hole and the redshift of the light which comes 
^v^j ' from the marginally trapped tube surface. 
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The study of black holes in stationary and asymptotically flat spacetimes has led to many remarkable 
^j 1 , insights. But, as we know, the our universe is not stationary and is in fact undergoing cosmological 

q-i expansion. Let us use the term cosmological black hole for any solution of Einstein equations representing 

a collapsing overdense region in a cosmological background-leading to an infinite density at its center 
(3JT)' [l|. The first attempt in this direction is due to McVittie who introduced a spacetime metric that 

represents a point mass embedded in a Friedmann-Robertson- Walker (FRW) universe. There have been 
several other attempts to construct solutions of Einstein equations representing such a collapsing central 
mass. Gluing of a Schwarzschild manifold to an expanding FRW manifold is one of such attempts, made 
, first by Einstein and Straus [3j. 

Now, a widely used metric to describe the gravitational collapse of a spherically symmetric dust cloud 
is the so-called Tolman-Bondi-Lemaitre (LTB) metric [j]. These models have been extensively studied 
for the validity of the cosmic censorship conjecture [1, @|- It was pointed out in Q that the model admit 
■ cosmological black holes. 

Much of the literature on BH's focuses on stationary and asymptotically fiat situation It is 
therefore desirable to have black hole models embedded in cosmological environment to see if there 
may be considerable differences from the familiar Schwarzschild black hole, and if there are quasi-local 
characteristics of it in place of the global and teleological event horizon. The need for local definition 
of black holes and their horizons has led to concepts such as Hayward's trapping horizons isolated 
horizons [lS^ . Ashtekar and Krishnan's dynamical horizon (DH)[l3j], and Booth and Fairhurst's slowly 
evolving horizon [l4j |. 

There have beeb some previous studies of BH's in cosmological situations, for example to glue two 
different LTB metrics to study the structure formation out of an initial mass condensation or the 
formation of a galaxy with a central black hole |l5| and [161 ] . However, since the structure of the metric 
outside these mass condensations are fixed by hand to match to a specific galaxy or cluster feature, we 
are faced with the shortcomings of the cut and paste models. Harada et al, was also interested in the 
behavior of primordial black holes within cosmological models to probe the gravitational memory and 
back hole radiation in expanding universe [l7i j . 

Our main interest is to consider the black hole (BH) properties within the expanding universe model. 
We consider different properties of LTB metric as a cosmological BH in section II. These properties are 
generalized to perfect fluid in section III. Section IV is devoted to cosmological constant effect on spherical 
cosmological BH. Throughout the paper we assume 8nG = c = 1. 
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II. LTB METRIC 

The LTB Q metric is a spherically symmetric nonstatic solution of the Einstein equations with a dust 
source. 

The LTB metric may be written in synchronous coordinates as 

ds 2 = -dt 2 + R ' dr 2 + R(t, r) W. (1) 
1 + f{r) 

and represents a pressure-less perfect fluid satisfying 

Here dot and prime denote partial derivatives with respect to the parameters t and r respectively. The 
angular distance R, depending on the value of /, is given by 

M 

R = -j(l-cos(r)(r,t))), 
f-fW 2 

v - s in( V )= { -^(t-t n (r)), (3) 



for / < 0, and 

R = ( 

2 

for / = 0, and 



R=&M)i(t -t n )i, (4) 



M 

r= (cosh(77(r,t)) - f), 

J3/2 

sinh(7?) -77= — -{t-t n (r)), (5) 

for / > 0. 

The metric is covariant under the rescaling r — > r(r). Therefore, one can fix one of the three free functions 
of the metric, i.e. t n (r), /(r), and M(r). One can shows that M(r) represents the mass accumulation in 
a 2-sphere with radius r, more precisely, the Misner-Sharp mass 0. 

There are two generic singularities of this metric: the big bang and big crunch singularity (shell focusing 
singularity) at R(t,r) = 0, and the shell crossing one at R'(t,r) — 0. However, if and is finite 

at R = then there is no shell focusing singularity. Similarity, if ^ is finite at R' = then there is no 
shell crossing singularity. In addition, to get rid of the unnecessary complications of the shell focusing 
singularity, corresponding to a non-simultaneous big bang singularity, we will assume t n (r) — 0. This will 
enable us to concentrate on the behavior of the collapse of an overdense region in an expanding universe 
without interfering with the complexity of the inherent bang singularity of the metric [7J. 
Before discussing about the horizon properties for LTB metric, we needs to know short summary of 
dynamical horizon. A marginally trapped tube T (MTT) is a hypersurface in a 4-dimensional spacetime 
that is foliated by two-surfaces S, called marginally trapped surfaces, such that #( n )|T < an d 6(i)\t = 0. 
MTTs have no restriction on their signature, which is allowed to vary over the hypersurface. This is a 
generalization of the familiar concept of the apparent horizon [l3| • If a MTT is everywhere spacelikc it is 
referred to as a dynamical horizon. If it is everywhere timelike it is called a timelike membrane (TLM). 
In case it is everywhere null and non-expanding then we have an isolated horizon. 

Proposition II. 1 There exists a marginally trapped tube for dust cosmological BH. 

As we know any dynamical BH has collapsing region (R < 0). If we calculate the expansion for ingoing 

and outgoing null geodesic then we get, 9^ oc (1 — /f + y )) #(n) * (~ 1 ~ vi+/ ^ < ®' -^ ence ' we see 
in collapsing region, the expansion for null outgoing geodesic changes its sign from negative to positive 
at R = 2M as we go to larger R and ingoing null geodesies expansion is negative everywhere. Therefore, 
the 3-manifold R = 2M is a marginally trapped tube (MTT). 

Now we prove that, there is MTT between singularity line (R = 0,77 = 2tt) and R = 0, rj = n (boundary 
between collapsing and expanding region). It can be seen from ([3]) that we have R > 2M on R = (it 
is sufficient to put / > — 1 and n = tt). We know that at the singularity (R=0) the Misner-Sharp mass 
in non zero M(r = r s ) — M s > R s — 0, so if we look at the R and 2M values at R plane, there is a i?o 
between R = and R = that i?o = 2M so we have the apparent horizon. This MTT can be spacelike, 
timelike or null surface. 
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A. LTB solutions as a cosmological BH 

For LTB solutions to be asymptotically FRW certain conditions have to be fulfilled. We first note 
that FRW spaces are special cases of LTB metrics: if R(r, t) is separated as R(r, t) = ra(t) we obviously 
get the homogeneous FRW solutions. For the vanishing bang time this corresponds to M(r) = cr 3 , 
f(r) — — r 2 , 0, r 2 . Therefore, to have an asymptotically open FRW LTB solution must have M(r) = cr 3 , 
f(r) = r\ 

Assumptions: The first assumption is that the Misner-Sharp mass of round spheres increases 
monotonically with respect to comoving coordinate r; M'(r) > 0. Our second assumption is that we 
have no shell crossing and R' > everywhere. The third assumption expresses the idea that want to 
look for conditions leading to an overdense region near the center r = 0, within an expanding universe 
with R > (at least) far from the center. However, overdensities in a region around the center require 
R < 0, corresponding to the collapse phase of the overdense region, which we may assume to start at 
a time t c > 0. We will say that a LTB solution represents a LTB cosmological BH if it satisfied these 
conditions. Finally, equations ([3]), and ([5]) it is easily seen that for the collapsing region one has 
to have f(r) < 0. In contrast, for the universe outside there must be f(r) > having expansion for 
late time (it can be f(r) < and we have expanding phase, but after some time this region change 
the phase form expanding to collapsing phase). Hence there must be at least one root for f(r) (see FigQJ. 

The collapse of the overdense region leads to two new conditions on the metric coefficients. First, we 
see from Q, that at any constant time shells corresponding to < rj(r) < ir are in an expanding phase 
and those corresponding to 7r < r)(r) < 2tt are in the collapsing phase. 

We will now look for LTB solutions fulfilling these assumptions. 



f(r) 



IV-FRW: k=-1 



Expanding region 



Collapsing region 



FIG. 1: Behaviors of the curvature function f(r) 



Proposition II. 2 There is a finite time for LTB cosmological BH at which ^\mtt becomes positii 
Proof: For LTB metric we have, 

R^^-^R + i-tRK-sL). 



(6) 



In apparent horizon, 4t|ah = b'-2M' 



According to before proposition. III. 1 1 , we have R\ah = 



+ f 7^ 0. We know form the above equation that 



f ■ M' f 

r' 2M'\ MTT = -jR + (-«*)( ^ - 3 y)- 



(7) 



The first term on right hand side can be negative or positive and the second term is +. If the first term be 
negative, it is easy to see from above equation that there will be a finite to that 4- J? < (— to-R)(^j- — 3y ) 
on MTT (because the second term is increasing according to t in collapsing region while the first term 



Singularity 

Apparent horizon 




FIG. 2: Causal structure for cosmological BH. 



remain finite collapsing region). Hence, %\ah > 0. 

As we see there is a finite time to a t which R' — 2M' > and so R' — M' > 0; therefor apparent horizon 
become space like (7| . This is a dynamical horizon which is good candidate for black hole boundary. 
We assume that the collapsing region is limited a finite region ro in the expanding region then apparent 
horizon is limited too. Consider following proposition from analysis, 

Proposition II. 3 If ^jp- > and r < ro, for any ei > and €2 > there will be to such that < t\ 

dt V t'2 ■ 

As t becomes larger and larger, the %\ah tends to zero. One can see that the r coordinate at event 
horizon has the same behavior of singularity line at late time [l8[ . 

Proposition II. 4 Singularity line slope int — r plane is increasing for cosmological dust collapse. 

We know that r is the comoving radius. We want to construct a collapsing spherical matter in expanding 
universe, so we expect that the singularity forms at the center, r = firstly and then any other shell fall 
into the singularity. Each shell which is more nearer to the center will fall into the singularity earlier. 
Therefore, the above expression int — r plane can be mathematically express as ^singularity = t' s ( r ) > 0. 
We should note that t s (r) and R = have the same behavior in t — r plane because they correspond to 
the 77 = 2tt and rj = ir respectively and their difference is only 2 factor. 

Now consider the first root of the f(r = ro) = which is outside collapsing region and is in the 
expanding phase (at any time). It has finite mass M which corresponds to 77 = 0. The singularity line 
t s (r) and R = cannot reach to these surface in any time at t — r plan. Hence, there is a r = ro (such 
that f(r = ro) = 0) that the t s (r) and R = curves become asymptote to r = ro in t — r plan at late time. 



Although the space time is not asymptotically flat, for the class of LTB BHs under consideration we 
can still introduce the notion at an event horizon as "the very last ray to reach future singularity" or 
"the light ray that divides those observer who cannot escape the future singularity from thus who do" 
then we have the following result. 

Proposition II. 5 Every trapped surface (T) for cosmological BH is located inside the event horizon. 

Our matter is dust and satisfies energy condition. Suppose that there is a point p that are located 
at outside the event horizon. Therefor, it can send a ray to point q infinity (not singularity). Hence, 
according to theorem 9.3.11 at [§] there is null geodesic 7 from p to q is orthogonal to T has no conjugate 
point between T and q. However, this is impossible, because according to theorem 9.3.6, 7 must have a 
conjugate point within affine parameter 2/|#o| fromp, where #0 < is the expansion at p of the orthogonal 
null geodesic congruence from T to which 7 belongs. 

Geometry of the dynamical horizon H is represented by The unit normal to H by r a ; g a bT a T b = — 1 . 
The unit space-like vector orthogonal to S and tangent to H is denoted by r a . Finally, we will fix the 
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rescaling freedom in the choice of null normals via l a = r a + r a and n a = r a — r a . As is usual in general 
relativity, the notion of energy is tied to a choice of a vector field. The definition of a dynamical horizon 
provides a preferred direction field; that along l a . To fix the proportionality factor, or the lapse N, let us 
first introduce the area radius R, a function which is constant on each S and satisfies as — AnR 2 . Since 
we already know that area is monotonically increasing, R is a good coordinate on H . Now, the 3- volume 
d 3 V on H can be decomposed as d 3 V = \dR\~ 1 dRd 2 V where d denotes the gradient on H. Therefore, 
as we will see, our calculations will simplify if we choose Nr = \dR\. 

Fix two cross sections Si and S2 of H and denote by AH the portion of H they bound. We are 
interested in calculating the flux of energy associated with = Nf>£ a across AH. Denote the flux of 

matter energy across AH by T^J ttcr : 

■p(. R ) — / 71 ,~atb j3ta _ 
•^matter - — / 1 ab~ ^(R) a v — 



matter — / J aii' 5(m ( 
J AH 

^(M(r 2 ) - M(n)). (8) 

LTB BHs, dynamical horizon behave like ([5J at late time, matter flux become zero and the horizon 
radius become fixed at radius R — 2M(ro) where M(tq) is Misner-Sharp mass in the horizon. For more 
practical case we can calculate the matter accretion according time along the dynamical horizon by 



dM(r) _ dM (r) dr 

-^r lDH -^r-d~t lDH - 



(9) 



To show that the dynamical horizon becomes isolated horizon, we have to show that the dynamical 
horizon become Slowly evolving horizon at late time. It can be seen from Q that 



dt, R' — 2M' R' - 2M' 



dr^ AH -R l + f 



(10) 



We know R' > has regular behavior (no shell crossing) and f(r) is finite along the MTT. From 
Proposition lIL3l we R' — 2M'\ah becomes large and R' — M'\ah do too. We follow the Kavanagh and 
Booth [19| calculations to show that our dynamical horizon becomes a slowly evolving horizon at late 
time. In defining a slowly evolving horizon, it is convenient to further restrict the scaling of the null 
vectors. To do this we label the foliating MTSs with a parameter v and choose the scaling and an 
evolution parameter C so that 

V = V 1 - cn^ (11) 

is tangent to DH and 

£ v v = l. (12) 
According to 19], if we calculate c and e along the DH for LTB metric, we have 

= 2^/1+7^^1^ (13) 



e 2 = 8v/iT7-^-^u ff (w) 

if c, e << 1, then slowly evolving horizons conditions become satisfied. As we see in above , R! — M'\ah 
becomes large at large t, therefor our DH becomes slowly evolving horizons at late time. We can see 
form the Einstein equation that the density become small around the apparent horizon (R and M 1 are 
finite because we have no singularity at there) and we will have a void around the apparent horizon at 
late time. 



B. Redshift 



If an emitter send a light ray to an observer with null vector k^, the relative light redshift that is 
calculated by observer with 4-velocity is, 



G 



From [2(|, we can can calculate the redshift for observer who sit at r = const with u M = (1,0,0,0) as 
below, 

fc* = c exp ^- y ^j_ rfr^ = 

•■*/**7(a-£+&*)- 

As we saw in the last section, i?' becomes large for late time in apparent horizon and {k^u^),, will have 
larger value as t become larger. Therefor, the light that come out from this surface to external observer 
will be fainted and has the infinite redshift relative to distance observer at large time. 

III. COSMOLOGICAL BH WITH PERFECT FLUID 

In last section we limited ourself to a dust source and studied at BH formation in expanding background. 
In this section we want to probe if the same properties of the dust cosmological BH collapse are hold for 
general perfect fluid collapse. We assume perfect fluid with an equation of state of the form p — wp. If 
w = the fluid becomes dust. We now assume that w ^ and has finite value. Take a collapsing ideal 
fluid within a compact spherically symmetric spacetime region described by the following metric in the 
comoving coordinates (t, r, 9, ip): 

dg 2 = _ e 2v(Ur) dt 2 + e 2^t,r) dr 2 + R ^ r) 2^2 

Assuming the energy momentum tensor for the perfect fluid in the form 

Tt = -p(t,r), T;=p r (t,r), T e e = 

T£=po(t,r)=wp(t,r), (18) 

with the week energy condition 

P > p + p r > p + pe > 0, (19) 
where w is constant. Einstein equations give, 

2M' 2M 



where 



and M is defined by 



R 

The function M can also be written as 



R*R" Pr = -lPR> m 



v , _ 2{p e -pr) R' P' r , 21 s 

P+Pr R P + Pr ' 



2R' + R'^+r¥-=0, (22) 



G = e- 2 ^{R') 2 , H = e~ 2 »(R) 2 , (23) 



2M 

G _ R = l_t^L. (24) 



1 r R 

M = - pR 2 dR, (25) 
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or 



where 



and 



tfV = Ane^R'dr, (27) 



#-=«"T- (28) 
dr dt y ' 

The last form of the function M indicates that when considered as energy, it includes contribution from 
the kinetic energy and the gravitational potential energy. M is called the Misner-Sharp energy. 



Hayward showed that in the Newtonian limit of a perfect fluid, M yields the Newtonian mass to 
leading order and the Newtonian kinetic and potential energy to the next order. In vacuum, M reduces 
to the Schwarzschild energy. At null and spatial infinity, M reduces to the Bondi-Sachs and Arnowitt- 
Deser-Misner energies respectively [ll|. Similar to dust case the flux of matter energy across AH by 
t-(«) . 

matter" 

-p( R ) — / T 1 py _ 

^matter - — / 1 ab~ <Z(R) a v — 
J AH 

i(M(t 2 ,r 2 )-M(t 1 ,r 1 )). (29) 



We assume M'(t, r) > and M(t, r) > in collapsing region which came form positive sign of density 
and pressure. 

Existence of marginally trapped tube for cosmological BH: As we know any dynamical BH has collapsing 
region (R < 0). If we calculate the expansion for ingoing and outgoing null geodesic then we get, 

1 2M ^ I 2M ^ 

6(i) oc (1 — y 1 H — q q — ), Q(n) (—1 — yl H — a g — ) < 0. Hence, we see in collapsing region, the 
expansion for null outgoing geodesic changes its sign from negative to positive at R = 2M as we go to 
larger R and ingoing null geodesies expansion is negative everywhere. Therefore, 3-manifold R = 2M is 
marginally trapped tube. 

For proof that there exist MTT between the singularity line and the boundary between ex- 
pansion region and collapsing region, we know that at R = so H = and from we have 
G > then R > 2m\ R _ , On the other hand in singularity 2M(t,r s ) > R = 0. If we look at the 
R and 2M values at R plane, there is a Rq between R — and R = that Rq = 2M so we have the MTT. 

Similarly, we assume that M' > (which means that Misner-Sharp mass is increasing ), we have no 
shell crossing e 2 ^' t,r ' > 0(or equivalently for finite G R 1 > 0), M > in collapsing region and M < 
in expanding region. In the region between collapsing and expanding matter; R = and then M = 
because from (|20p we have no singularity in pressure. These condition are compatible with p and p 
positive sign. 



Proposition III.l the Singularity line slope int — r plane increases for general spherical fluid collapse. 

We want to show if we assume that singularity line has other shape for example part A and part B ( 
see FigJ3]), then we get contradictions with our assumptions. Consider part A, if these case exist then we 
have R{r\) > Rfa) = which is in contradiction with R' > and Consider part B, if these case exist 
then we have R(t\) > Rfo) = which is in contradiction with R < 0. Therefor, t' s > comes from the 
above assumptions for cosmological BH collapse. 

Cosmological expansion assumption force singularity to be limited in the part of of the space time 
around the center, and the other part expand with cosmological expansion. Hence, there is an event 
horizon that separates causal geodesies that fall in to the singularity from other geodesies. This event 
horizon is asymptotic to the singularity at late time. Similar to the dust case, the trapped region is inside 
event horizon see Proposition lII.5l 



Prnnnsitinn TTT.2 The n.nnn.rent hnriznn. slnnp. in, t — r nln.n.e hp.cnm.ps infinite n.t late tim.p. 
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FIG. 3: Impossible shape for the singularity line in BH collapse. 



We know that the singularity slop is strictly positive and singularity is limited in a finite r, therefor its 
slop become infinite (see Proposition lII.3| ). The apparent horizon is among the singularity line and the 
event horizon, so as even horizon tend to the singularity then apparent horizon tend to singularity and 
it becomes asymptotic to singularity at late time and its slop becomes infinite. 
We want to show that MTT becomes slowly evolving horizons large time. We know the expression 

dt. R' — 2M' 

is large at from above proposition late time. We know R' > has regular behavior (no shell crossing for 
finite G), R < and M > are finite along the apparent horizon. Therefor, we can infer 2M — R is 
finite. Since the density is finite around the horizon, so R 1 — 2M' = R' — 2R'pR 2 should be large and 
then R'\ah 3> becomes large. 

If we want than MTT becomes space like, we should have 

J^L = -i< [t *U<i. (31) 

dr l nul1 \ 1 J" / 

Because of M\ah > and R\ah < 0, the above condition will be satisfied and MTT becomes dynamical 
horizon at large time. 

Similar to LTB case, it can be shown that the dynamical horizon satisfying slowly evolving horizons at 
large time. According to [llj], if we calculate c and e along the DH for perfect fluid metric, we have 

M'R + wM'R , 

C = 2 ; : ; \AH (32) 

M'R - wM'R - RR' 

£2 = g M>R + vM>R 

M'R — wM'R — RR' 

which we put M — —wM'jp from the perfect fluid assumption. If c, e << 1, then slowly evolving 
horizons conditions become satisfied. As we calculate in above , R'\ah become large at large t, therefor 
the dynamical horizon becomes slowly evolving horizons at large time. Similar to LTB, the density 
becomes small around the {R and M' are finite around the MTT) and void will be formed around the 
dynamical horizon at late time. 



IV. COSMOLOGICAL BH WITH COSMOLOGICAL CONSTANT 



One of the characteristic properties of a black hole is its horizon. One can show that cosmological 

constant ran Via.vp an pffprt on thp na.tnrp of thp sincmlaritv and final fatp of trip VilacV Violp l2l! Rut wp 
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want to show that local behavior of the gravitational collapse are not affected by cosmological constant 
as a dark energy candidate. If we calculate the expansion for ingoing and outgoing null geodesic then we 

/ AR 2 | 2M | f / Aft 2 | 2M | j: 

get, 8^ oc (1 — — — 3 ^ 1+ J )> ( — 1 — - — 3 yi+f ^ ^ ®' ^° tne a PP aren * horizon surface change 

from 2M = R to 2M = R — ^j- surface with the same properties. 

For general perfect fluid we have the same change in apparent horizon surface from expansion of light 

/ Ai? 2 | 2M ^ I AR 2 | 2AI ^ 

ray, 8^ cx (1 — y 1 H - — ^ ), #(„) oc (—1 — y H 3 — ^ ) < (The apparent horizon surface 

2M = R-^). 

^ttcr := J AH (T ab - Jt_g ab )T*e (R) d 3 V. (34) 
It can be seen from fijl that the area law becomes, 



dE* =^-{l- KR 2 )d(AnR 2 ) (35) 



which kr — j£ is effective surface gravity. We can integrate from the above equation 

If we compare a cosmological BH with a radius i?o in de Sitter background (A > 0) with non de Sitter 
background, we will get that de Sitter background black hole will have less Misner-Sharp mass M in 
black hole with radius Rq. 

The local property of the black hole evolution will not change with adding the cosmological constant 
term because if we put cosmological constant realistic value in the Einstein equation with A^; 10~ 35 s~ 2 
and mass M = 1O 12 M , we get 

2M A , A 9 2M 
G-H = l-— + -R 2 , -i? 2 «— , (37) 

jR 2 < 10~ 20 . Hence, cosmological constant don't affect on local black hole physics. 



V. DISCUSSION 



Considering the back hole as a dynamical object in expanding background shows new properties that 
a schwarzschild back hole is not able to show. In particular, the mass and the flux of matter can have 
different properties in expanding background. Recent analysis has shown that these dynamical BHs have 
flexibility to explore some basic questions such as cosmic censorship conjecture, quasi-local definition 
masses and BH thermodynamic in non asymptotic flat back ground [H Il3j. Thus, the fist attempt 
is making a good model for these BH (cosmological BH) with reasonable kind of matter Q. In this 
paper we assumed some acceptable properties for these BH and consider its singularity and its horizons 
behavior along the time. We found that the singularity line shows a special behavior and MTT becomes 
dynamical horizon as a candidate for BH boundary in non asymptotic flat background. Furthermore, 
in terms of global behavior in time, Misner-Sharp mass has this ability to describe the mass of the 
cosmological BH and its time evolution across the dynamical horizon becomes the fluid accretion. In 
contrast the schwarzschild BH, the light with infinite redshift cannot come from the apparent horizon 
to distance observer. We showed that for the dust cosmological collapse, the apparent horizon becomes 
infinite redshift surface only at late time for distance observer. Finally, we found that we can neglect the 
cosmological constant (as a dark energy candidate) for local behavior of the BH in time. 
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